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Abstract 

The smallness of the observed neutrino masses might have a radiative ori- 
gin. Here we revisit a specific two-loop model of neutrino mass, independently 
proposed by Babu and Zee. We point out that current constraints from neu- 
trino data can be used to derive strict lower limits on the branching ratio of 
flavour changing charged lepton decays, such as ;U — > 67. Non-observation of 
Br(^ — > 67) at the level of 10~^^ would rule out singly charged scalar masses 
smaller than 590 GeV (5.04 TeV) in case of normal (inverse) neutrino mass 
hierarchy. Conversely, decay branching ratios of the non-standard scalars of 
the model can be fixed by the measured neutrino angles (and mass scale). 
Thus, if the scalars of the model are light enough to be produced at the LHC 
or ILC, measuring their decay properties would serve as a direct test of the 
model as the origin of neutrino masses. 



1 Introduction 



During the past few years neutrino oscillation experiments have firmly established 
that neutrinos have non-zero masses and mixing angles among the different genera- 
tions pp. While for the absolute scale of neutrino mass only upper limits of the order 
nil, ~ C>{2 eV) exist j2], two neutrino mass squared differences and two neutrino an- 
gles are by now known quite precisely jH] . These are the atmospheric neutrino mass, 
Aml,^ = (2.0 - 3.2) [IQ-^ eV^], and angle, sin^ ^Atm = (0.34 - 0.68), as well as the 
solar neutrino mass Am| = (7.1-8.9) [10"^ eV^], and angle, sin^ 6^0 = (0.24-0.40), 
all numbers at 3 a c.l. For the remaining neutrino angle, the so-called Chooz |1] or 
reactor neutrino angle dpt, a global fit to all neutrino data jB] currently gives a limit 
of sin^ Or < 0.04 @3a c.l. 

From a theoretical perspective, there exist several options to explain the small- 
ness of the observed neutrino masses. Perhaps the simplest - but certainly the most 
popular - possibility is the seesaw mechanism OIHl- Many variants of the seesaw 
exist, see for example the recent review [7j. However, most realizations of the seesaw 
make use of a large scale, typically the Grand Unification Scale, to suppress neutrino 
masses and are, therefore, only indirectly testable. 

On the other hand, many neutrino mass models exist, in which the scale of 
lepton number violation can be as low as the electro-weak scale or lower. Examples 
are supersymmetric models with violation of R-parity jH| E]; models with Higgs 
triplets or a combination of both [2] , leptoquarks ^2] or radiative models, both 
with neutrino masses at 1-loop [131 E] or at 2-loop [121 US Hn] order. Radiative 
mechanisms might be considered especially appealing, since they generate small 
neutrino masses automatically, essentially due to loop factors. 

In this paper we will concentrate on a model of neutrino masses, proposed inde- 
pendently by Zee and Babu [inj, in which neutrino masses arise only at 2-loop 
order. The model introduces two new charged scalars, and k^~^, both singlets 
under SU{2)l, which couple only to leptons. One can easily estimate, see fig. (HJ 
and the discussion in the next section, that neutrino masses in this setup are of order 
nil, ~ {ph)/{16TT'^)^{m'j^/'ms), i.e. (9 ~ 1 eV for couplings / and h of order C(l) 
and scalar mass parameters, mg, of order (9(100) GeV. Given that current neutrino 
data requires at least one neutrino to have a mass of order (9(0.05) eV, one expects 
that the new scalars should have masses in the range (9(0.1 — 1) TeV. The model is 
therefore potentially testable at near-future accelerators, such as the LHC or ILC. 

Babu and Macesanu J7| recently re-analyzed this model in light of solar and 
atmospheric neutrino oscillation data. They identified the regions in parameter 
space, in which the model can explain the experimental neutrino data and tabulated 
in some detail constraints on the model parameters, which can be derived from the 
non-observation of various lepton flavour violating decay processes. Here, we extend 
upon the results presented in [T7| by pointing out that (a) current neutrino data 
can be used to derive absolute lower limits on the branching ratios of the processes 
^f3l- Especially important in view of future experimental sensitivities |18| is 
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Figure 1: Feynman diagram for the 2-loop Majorana neutrino masses in the model 

of [laiin]. 

that Br{fi 67) > 10^^^ is guaranteed for charged scalar masses smaller than 
590 GeV (5.04 TeV) in case of normal (inverse) neutrino mass hierarchy. And (b) 
decay branching ratios of the non-standard scalars of the model can be fixed by the 
measured neutrino angles (and mass scale). Thus, if the scalars of the model are 
light enough to be produced at the LHC or ILC, measuring their decay properties 
would serve as a direct test of the model as the origin of neutrino masses. 

The rest of this paper is organized as follows. In the next section, we discuss 
the Lagrangian of the model, as well as its parameters in light of current oscillation 
data. In this part we will make extensive use of the results of ^7j. In section 
3, we calculate flavour violating charged lepton decays, la — > h^Jd and l/^'j, 
discussing their connection with neutrino physics in some detail. Then, we consider 
the decays of the new scalars at future colliders, presenting ranges for various decay 
branching ratios as predicted by current neutrino data. We then close with a short 
discussion. 



2 Neutrino masses at 2-loop 

As mentioned above, the model we consider [TS1|^ is a simple extension of the stan- 
dard model, containing two new scalars, h'^ and k^~^, both singlets under SU{2)l. 
Their coupling to standard model leptons is given by 

£ = U{LllCL^p,)e,,h^ + K^{el^Cepn)k^^ + h.c. (1) 

Here, are the standard model (left-handed) lepton doublets, cr the charged 
lepton singlets, a,j3 are generation indices and e^j is the completely antisymmetric 
tensor. Note that / is antisymmetric, while h' is symmetric. Assigning L = 2 to h~ 
and k~^^, eq. (P) conserves lepton number. Lepton number violation in the model 
resides only in the following term in the scalar potential 

C = -fxh^h+k— + h.c. (2) 

Here, is a parameter with dimension of mass, its value is not predicted by the 
model. However, vacuum stability arguments can be used to derive an upper bound 
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for this parameter ^Tj. For ruh ~ this bound reads 

< (67r2)VV. (3) 

The setup of eq. ((T)) and eq. Q generates Majorana neutrino masses via the two- 
loop diagram shown in fig. (P). The resulting neutrino mass matrix can be expressed 
as ^ 

with summation over x, y implied. The parameters defined as tUxy = 

mxhxyfny, with rrix the mass of the charged lepton Ix- Following ^Tj we have rewrit- 
ten haa = h'^a ^a/s = X(r) finally is a dimensionless two-loop integral 
given by ^ 

^r) = -fdx f'^ dy ^ log y^^pl. (5) 

Jo Jo X + [r — l)y + y^ x + ry 

For non-zero values of r, X(r) can be solved only numerically. We note that for the 
range of interest, say 10~^ < r < 10^, X(r) varies quite smoothly between (roughly) 
3 < X(r) < 0.2. 

Eq-® generates only two non-zero neutrino masses. This can easily be seen from 
its index structure: Det(A^'') = Get^fax^xyfy/s) = Det(/a/3) = 0. The model there- 
fore can not generate a degenerate neutrino spectrum. One can find the eigenvector 
for the massless state, it is proportional to 

v^ = mh-e,e') (6) 

where A/" = (1 + + e'^)~^/^ is a normalization factor. Here we have introduced 

J J IJ,T 

With Ai'^.VQ = one can express the parameters e and e' also in terms of the entries 
of the neutrino mass matrix. A straightforward calculation yields 

mumss - mi3m23 



^12^23 - ^13^22 



(8) 



Interestingly, eq. (jH)) can be rewritten directly as a function of the measured neutrino 
angles. For normal hierarchy, i.e. m^^^^ ~ (0,m,M), one obtains ^ 

e = tan^l2^^^^^ + tan^l3sin^23e"^^ (9) 

COSt^l3 

e = tan6'i2— — ^ — tan6'i3 cos^23e~*''. 
cos o'ls 



^We correct a minor misprint in eq.(7) of jl7| . 



^We use the notation m ~ y Atoq and M ~ y^Am^^, as well as m^^ ~ for inverse hierarchy. 
This has the advantage that 9i2 = 6q, 623 — Oauh and 613 = Oji for both hierarchies. 
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Note, that eq. Q does not depend on neutrino masses, and that current data on 
neutrino angles require both e and e' to be non-zero. On the other hand, in the case 
of inverse hierarchy, j 3 ~ (M, ±M + m, 0), eq. (jH)) leads to 



— cot sin 6*236' 

cot ^13 cos ^2 



(10) 



Again, both e and e' have to be different from zero. Note that 6 in eq. and 
(ll(J|) is a CP-violating phase, which reduces to a CP-sign 5 = 0, tt in case of real 
parameters. 

With the equations outlined above, we are now in a position to give an estimate 
of the typical size of neutrino masses in the model. For an analytical understanding, 
the following approximation is quite helpful. Since me -C m^,m^, cj, 
are expected to be much smaller than the other 00^/3- Then, in the limit = u, 
ojer = 0, eq. (jH) reduces to 



Uef, and 



I 



'111 



where 



c 



.ml 



(167r^)2 m 



(12) 



From eq. (fTTj) it is easy to estimate the typical ranges of parameters, for which 
the model can explain current neutrino data. In case of normal hierarchy, a large 
atmospheric angle requires uj 



TT- Thus, we find the constraint 



m 



1 



(13) 



^ requires e ~ e' ~ 1/2, see 
~ uJt-ti although with a different 

M 



— , see also eq. 



On the other hand, a solar angle of order tan^0 
eq. ©• Inverse hierarchy still requires ^ UJ^^ 
relative sign, while e and e' have to be much larger, i.e. e 

m- 

What is the maximal neutrino mass the model can generate? Using eqs Q 
and (jl3p . this upper limit can be estimated choosing h^^ maximal. Motivated by 
perturbativity, we choose h^^ = 1. ^ Then, > 800 GeV is required (see the next 
section), and with rrih = 100 GeV, X(r) ~ 0.3 results. Putting finally f^^ = 0.03 we 
arrive at rrij^g ~ 0.05 eV. Since all other parameters in this estimate have been put to 
extreme values, Z^,- > 0.03 will be required in general. Obviously, even considering 
only neutrino data, the parameters of the model are already severely constrained. 

■^One could also choose /i^^ = y/in . However, as pointed out in ^7\, even /i^^ = 1 at the weak 
scale will result in non-perturbative values of /i^^ at scales just one order of magnitude larger. 
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3 Flavour violating charged lepton decays 



Due to the flavour off-diagonal couplings of the k^~^ and scalars to SM leptons, 
the model has sizeable non-standard flavour violating charged lepton decays. An 
extensive list of constraints on model parameters, derived from the observed upper 
limits of these decays, can be found in [T2j . Here we will discuss decays of the type 
^/37 and their connection with neutrino physics. As the experimentally most 
interesting case we concentrate on /i — >■ ej. A short comment on r decays is given 
at the end of this section. 

Consider the partial decay width of la — > l/s'j induced by the scalar loop 
shown in fig. (0). In the limit of <^ it is given by 



T{L - hi) = 2am^(-5)^p^ . (14) 

We will be interested in deriving a lower bound on the numerical value of eq. ()14|) in 
the following. Note, that although there is a graph similar to the one shown in fig. 
(j2I) with a k'^^ in the intermediate state, there is no interference between the two 
contributions (in the limit where the smaller lepton mass is put to zero). Thus, in 
deriving the lowest possible value of Br(/i — > 67) we will put the contribution from 
fc"*""*" to zero. Any finite contribution from the doubly charged scalars would lead to 
stronger bounds on vrth than the numbers quoted below. 
Using eqs ((Tj), (fTT|) and (|T^ we can rewrite eq. (fT^ as 

^ ^ 18v^r^/i2^J(r)2m2 ^ ^ 

^ 4.5-10-^°( w^n^AVl00GeV^2 

V/i2^j(^)2no.05eV^ ^ rriH > ^ ' 

With e non-zero, constrained by eq. or eq. (fTIHl in case of normal or inverse 
hierarchy, Br(/x —>■ 67) has to be non-zero as well. Its smallest numerical value is 
found for the largest possible value of /i^^ and X(r). 

In order to calculate X(r) we need to fix r consistent with all experimental 
constraints. This is done in the following way. The decay width — > IblJd induced 
by virtual exchange of k^^, see fig. (j21), is, in the limit rnhjincmii <^ m^, 

r(/, ^ = 1 Jli- \ (17) 

The most relevant constraint for the current discussion is derived from the upper 
bound on r — i> 3/i decay, which yields. 



\hnrh 



ml 



< GeV-\ (18) 
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Figure 2: Example diagrams for flavour changing charged lepton decays in the 
model. In addition to the diagrams shown, there are also box graphs involving 
contributing to la — ^ hlJ-d-, as well as graphs with k^^ , similar to the one shown, 
contributing to la hi- 



Br[ii —J- 67) 5r(/i — > 67) 




ruh [GeV] mu [GeV] 



Figure 3: Conservative lower limit on the branching ratio Br(/i — > 67) as a function 
of the charged scalar mass mh for normal hierarchy. The three lines are for the 
current solar angle sin^ 612 best fit value (full line) and 3 a lower (dashed line) and 
upper (dot-dashed line) bounds. To the left 5 = 0, to the right 5 = n. Other 
parameters fixed at sin^ ^23 = 0.5, sin^ ^13 = 0.040 and /\m\^^ = 2.0 ■ 10~^ eV^. 
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Figure 4: Dependence of the lower limit on Br(yU — > 67) for normal hierarchy on 
neutrino angles, for /i^^ = 1,6 = vr, m/^ = lO'^GeV. Left plot: (sin^6'23 = 0.68, 
sin^^is = 0) dashed line, (sin^^ss = 0.68, sin^^ig = 0.040) full line, (sin^ ^23 = 0.34, 
sin^^is = 0) dash-dotted line, (sin^ ^23 = 0.34, sin^ ^13 = 0.040) dotted line. Right 
plot: (sin2 6'i2 = 0.40, sin^ 6I13 = 0) dash-dotted line, (sin^ 6*12 = 0.40, sin^ 6113 = 
0.040) dotted line, (sin2^i2 = 0.24, sin^ ^13 = 0) dashed line, (sin^ ^12 = 0.24, 
sin^^is = 0.040) full line. 



For /i^r(^) = = 1, this bound implies nik > 770 GeV. For any fixed value of 
h^^, we can therefore fix the minimum value of r, i.e. the maximum allowed value 
of X(r), which in turn fixes the lower bound on Br(/i — 67). 

Fig. shows the resulting lower limit on Br(/x — 67) as a function of the 
charged scalar mass ruh for the case of normal hierarchy. In this plot, we have 
assumed that the parameters /i, /i^^ (and Am^^^^) take their maximal (minimal) 
allowed values, thus we consider this limit conservative. We would like to stress 
again, that any non-zero contributions to the decay /i — > 67 from k~^~^ can only 
increase Br(/i — > 67). 

Fig. (jH) and (0) show the dependence of the limit on Br(yU — > ej) on the three 
neutrino angles. Both plots are for the case of normal hierarchy. Larger values of 
^12 (^'23) result in larger (smaller) upper bounds. Smaller ranges of these parameters 
obviously lead to tighter predictions. For 6*13, below approximately sin^6'i3 < 0.01 
the dependence of Br(/i 67) is rather weak. 

Fig. (jni) shows the calculated lower limit on Br(/i 67) for the case of inverted 
hierarchy, both, versus the reactor angle and versus rrih. Due to the fact that e must 
be larger than e ~ sin6'23/ tan6'i3, the expected values for Br(/i — > 67) turn out to 
be much bigger than for the case of normal hierarchy. Even Br(/i — > 67) < 10^^^ 
requires already TeV-ish masses for m/j. 

The most conservative limits for ruh are always found for 5 = tt, sin^ 612 = 
(sin^^o)^'^^", sin2^23 = (sin^ ^Atm)*'^"'' and sin^ ^13 = (sin^ ^r)^^'*^. For the current 
bound of Br{fi ^ 67) < 1.2 x 10"", we find mh > 160 GeV {mh = 825 GeV) for 



7 



Br{ii 67) X 10 

4.0 



1.0 
0.6 

0.2 

10"^ 10"^ 

sin^ 

Figure 5: Dependence of the lower limit on Br(/i — > ej) for normal hierarchy on 
the reactor angle, for /i^^ = 1,6 = Ti.,mh = lO^GeV. Other parameters are chosen 
as (sin^^ss = 0.68, sin^ ^12 = 0.40) dashed line, (sin^ ^23 = 0.68, sin^ ^12 = 0.24) 
full line, (sin2^23 = 0.34, sin^ ^12 = 0.40) dash-dotted line and (sin^^sa = 0.34, 
sin2 6'i2 = 0.24) dotted line. 

normal (inverse) hierarchy. Future experiments expect to lower this limit to 
Br{fi — > 67) < 10"^'^, resulting in nih > 590 GeV (nih = 5040 GeV). Given these 
numbers, one expects that the MEG experiment ^18] will see the first evidence for 
fi —>■ ej in the near future, if the model discussed here indeed is the origin of neutrino 
masses. 

Finally, we would like to mention that the decays r ^ /i7 and r — 67 can 
be constrained in a similar way. However, the resulting lower limits, also of or- 
der 0(1O~^^), are far below the near-future experimental sensitivities and thus less 
interesting. 

4 Accelerator tests of the model 

In this section we will briefly discuss some possible accelerator signals of the model. 
With the couplings of h'^ and k^~^ tightly constrained by neutrino physics and 
flavour violating lepton decays, it turns out that various decay branching ratios can 
be predicted. Observing the corresponding final states could serve as a definite test 
of the model as the origin of neutrino masses. 

In fTTJ it has been estimated that at the LHC discovery of k^^ might be possible 
up to masses of < 1 TeV approximately. In the following we will therefore always 
assume that < 1 TeV and, in addition, rrih < 0.5 TeV. Given the discussion of 
the previous section, this range of masses implies that /i — ^ 67 should be seen at the 
MEG experiment. 
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Figure 6: Lower limit on Br(/i — > 67) for inverted hierarchy, to the left: versus the 
reactor angle; to the right: versus rrih. Parameter choices as before. The three lines 
are for the current sin^ 623 best fit value (full line) and 3 a upper (dot-dashed line) 
and lower (dashed line) bounds. 



The will decay to leptons with a partial decay width of, in the limit rria = 



1671 



(19) 



We can re-express eq. ()19p in terms of the parameters e and e' as 

l + e'2 



(20) 



^^(^'^^5'^^) = 2(l + e2 + e-)' 

It is therefore possible to directly "measure" or e'^ by calculating ratios of branch- 
ing ratio differences, such as the ones shown in fig. (|7j). Here, 



Brii'' = Br{h' 



Br{h- 



Br{h' 



(21) 



The plots on the left in fig. (|7j) show calculated ratios of branching ratios versus eq. 
Q, i.e. normal hierarchy, versus e (top) and e' (bottom). All points are obtained by 
numerically diagonalizing eq. (0} for random parameters and checking for consis- 
tency with all experimental constraints. However, since ^13 is unkown, eq. Q can 
be numerically calculated, but at the moment not experimentally determined. Thus, 
the plots on the right of the figure show the same ratios of branching ratios, but ver- 
sus (tan ^12 cos 6*12)^ and (tan ^12 sin Qyi)^ . The cut on sin^ 6'i3 of sin^ ^^13 < 2.5 x 10^^ 
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Figure 7: Ratios of decay branching ratios Br^^ , see text, versus /i (6*12) ^23) ^13) = 
tan ^12^ + tan ^13 sin 023 (top); and /2(0i2, ^23, ^13) = tan0i2f||| -tan^is cos023 
(bottom). In the right plots sin^ ^13 < 2.5 x 10^^ has been assumed. All points 
satisfy updated experimental neutrino data. 

in this plot is motivated by the expected sensitivity of the next generation of reactor 
experiments 120] • The width of the band of points in these plots indicates the 
uncertainty with which the corresponding ratios can be predicted. 

In case of normal (inverse) hierarchy, assuming best fit parameters for the neu- 
trino angles, eq. (j^n|) indicates that the branching ratios for e, /x and r final states 
of h+ decays should scale as 2/12 : 5/12 : 5/12 (1/2 : 1/4 : 1/4). Inserting the 
current 3 a ranges of the angles, following eqs. (jHl) and ()10|1 results in the following 
predicted ranges 

Br{h+ -> e ^ z/^) = [0.13, 0.22] ([0.48, 0.50]) (22) 
Br{h+ ^ fiJ2u/3) = [0.31,0.50] ([0.17,0.34]) 
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Br{h+ ^r^z/^) = [0.31.0.50] ([0.18,0.35]) 

/3 

for normal (inverse) hierarchy. The different predicted branching ratios for final 
states with electrons should make it nearly straightforward to distinguish normal 
and inverse hierarchy. Measuring any branching ratio outside the range given in eq. 
would rule out the model as possible origin of neutrino masses. 
The doubly charged scalar of the model decays either to two same-sign leptons 
or to two final states. The partial width to leptons is, for ma,mp = 0, 



r(fc++ - Lh) = ^hl, (23) 



whereas the decay width to two can be calculated to be 



r(fc++ ^ h^h^) = T^^/?(4) (24) 



Here, (i{x^) = y/T—Ax^ is a kinematical factor. 

The couplings h^is in eq.ipHj) are constrained by neutrino physics, see eq. lfT^ . 
and by lepton flavour violating decays of the type — > hh^d- For < 1 TeV the 
couplings hee, hen ^"^^ ^er are constrained to be smaller than 0.4, 4- 10~^ and 7- 10~^ 
[T7j . Thus, the leptonic final states of k~^~^ decays are mainly like-sign muon pairs 
(and possibly electrons). 

An interesting situation arises, if > 2mh. In this case, one can measure the 
lepton number violating parameter fi of eq.© by measuring the branching ratio of 
k^^ h^h^. Combining eq. and eq. we can write 

Br{k^^ h^h^) ^ ^ J"^}^ (25) 

Here, h^e <^ h^^ has been assumed. (For non-zero hee replace simply /i^^ hfj,f^ + hee 
in eq. (j2Hl)-) Plots of constant Br(A;++ h'^h'^) in the plane (m^, rrih) are shown in 
fig. Here, /i = frrih, with / = (Gvr^)^/^ has been used. 

Fig. (jSI) shows the resulting branching ratios for 2 values of /i^^, fixing in both 
cases the couplings /^^ such that the atmospheric neutrino mass is correctly repro- 
duced. For h^n < 0.2 the current limit on Br(/i —>■ 67) rules out all m/j < 0.5 TeV, 
thus this measurement is possible only for /i^^ > 0.2. Note that smaller values of n 
lead to smaller neutrino masses, thus upper bounds on the branching ratio for Br^'^ 
can be interpreted as upper limit on the neutrino mass in this model. 



5 Conclusion 

The observed smallness of neutrino masses could be understood if it has a radiative 
origin. In this paper, we have studied some phenomenological aspects of one well- 
known incarnation of this idea |15| ll6j. in which neutrino masses arise only at 2- loop 
order. 
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Figure 8: Lines of constant Br(A;++ h^h'^), assuming to the left /i^^ = 1: 
Er^'* = 0.1,0.2,0.3 and 0.4 for dotted, dash-dotted, full and dashed line. The 
vertical line corresponds to rrih = 208 GeV for which Br{n — > 67) = 1.2 x 10~^^ and 
horizontal line to = 743 GeV for which Br{T — 3fi) = 1.9 x 10~®, i.e. parameter 
combinations to the left/below this line are forbidden. Plot on the right assumes 
h^^ = 0.5. Lines are for Br^^ = 0.4, 0.5, 0.6 and 0.7, dotted, dash-dotted, full and 
dashed line. Again the shaded regions are excluded by Br(/i — 67) and Br(r — > 3/i). 



Given the observed neutrino masses and angles, it turns out that the parameters 
of this model are very tightly constrained already today and thus it is possible to 
make various predictions for the near future. Perhaps the phenomenologically most 
important one is, that one expects that the process /i ^ 67 has to be observed 
in the next round of experiments, i.e. Br(/i 67) > 10^^^ is guaranteed for 
singly charged scalar masses smaller than 590 GeV (5.04 TeV) for normal (inverse) 
hierarchical neutrino masses, and larger or even much larger branching ratios are 
expected in general. At least for the case of inverse hierarchy an upper limit on the 
decay — > 67 of this order would certainly remove most of the motivation to study 
this model. 

On the other hand, if fi —>■ ej is observed in the near future, it will be interesting 
to search for the charged scalars of the model at the LHC. As we have shown, in 
this case, several branching ratios of the decays of both, the singly and the doubly 
charged scalar are tightly fixed, mainly by data on neutrino angles. Observation 
of branching ratios outside the ranges discussed, would then definitely rule out the 
model possible explanation of neutrino masses. 
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